Current Reversals and Synchronization in Coupled Ratchets 
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Current reversal is an intriguing phenomenon that has been central to recent experimental and 
theoretical investigations of transport based on ratchet mechanism. By considering a system of 
two interacting ratchets, we demonstrate how the coupling can be used to control the reversals. In 
particular, we found that current reversal that exists in a single driven ratchet system can ultimately 
be eliminated in the presence of a second ratchet. Furthermore, we show that for larger values of 
the coupUng strength, the dynamics of the single driven ratchet completely dominates that of the 
coupled ratchets, thereby resulting in a fully synchronized state. 

PACS numbers: 05.45.Xt, 05.45.Ac, 05.40.Fb, 05.45.Pq, 05.60.Cd 
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Transport phenomena and, particularly, directed 
transport occur in many situations ranging from physi- 
cal systems to chemical and biological systems. Some re- 
cent research interest in transport problems is related to 
ratchet physics where unbiased, noise-induced transport 
occurs away from thermal equilibrium as a result of the 
action of Brownian motors [l|,|2|,0]. Brownian motors, es- 
pecially "ratchet" models, have been widely investigated 
partly due to the challenge to describe and control mech- 
anisms of fundamental biological processes at both the 
cell level (e.g. transport in ion channels) and body level 
(muscle operations) Another motivation is derived 
from recent advances in technology wherein devices for 
guiding tiny particles on nano/micro scales are sought; 
these include particle separation techniques, smoothing 
of atomic surfaces during electromigration and control of 
the motion of vortices in superconductors [H, [Hj] . Remark- 
ably, experimental realizations of some of these practical 
systems have been reported [6|, 0, HI • 

In this framework, two basic types of ratchet models 
have commonly been employed, namely: (i) the rock- 
ing or correlated ratchet, in which the particle is subject 
to an unbiased external force or additive noise, and (ii) 
and the flashing ratchet, where in the particle is periodi- 
cally kicked. The vast majority of these models are over- 
damped where the noise plays a vital role in the transport 
process. However, recent studies have shown that the role 
of noise can be replaced by deterministic chaos induced 
by the inertial term [9]. In such inertial ratchets, the 
issue of current reversal has been carefully investigated 
Besides, Hamiltonian ratchets have recently seen 
a breakthrough in the ratchet community [12]. 



In the description of the model transport in general, 
attention has mostly been paid to single particle ratch- 
ets. However many such systems coexist in numbers and 
do work in cooperation. For instance, molecular motors 
do not operate as a single particle but in groups - the 
most prominent example being the actin-myosin system 
in muscles [13]. For this reason, the relevance of many 
interacting particles and possible effects of collective be- 
havior, e.g. transport enhancement, current inversion, 
clustering and spontaneous current among others have 
been considered in the overdamped case ( UM a nd refer- 
ences therein) and the underdamped caseflBI, [l6l , [13] • 
Besides a few studies [l5|, [l6|, [l7|] on synchronization be- 
havior and possible connections to transport properties 
in the underdamped cases, a large number of studies in- 
vestigate transport properties based on the coupled over- 
damped case. 

In this paper, we study the synchronization dynamics 
of two coupled driven ratchets and show how the cou- 
pling can be used to control current reversals. We further 
show that the process of transition to full synchrony is 
preceded by several sudden changes in the current and 
particularly current reversals. 

The model we are interested in is made of two rocking 
ratchets [10], symmetrically perturbed by means of an 
elastic coupling. Its dynamics (in dimensionless form) 
can be described by 



hxj 



dV{xi,X2) 
dxi 



acos{ujt) (z = l,2), (1) 
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where the normalized time t is taken in units of the small 
resonant frequency uJo~^ of the system; a and uj represent 
the amplitude and the frequency of the driving, respec- 
tively; and b the damping parameter. Here V{xi^X2) 
is the perturbed two-dimensional ratchet potential given 
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FIG. 1: (Color online) Equipotential contours plot of 
y(xi,X2) with colors growing from blue (minima) to red 
(maxima): (a) No interaction, /c = 0, (b) weak coupling, 
k = 0.05, (c) moderate coupling, k — 0.15 and (d) strong 
coupling, k = 1.0. 



as: 

Vixi,X2) = 2C - ^ [<PiXi) + <j>iX2)] + X2f, 

(2) 

where ^(x^) = sin27r(xi — Xq) + |sin47r(x^ — Xq), i = 
1, 2 and k is the coupling strength which determines the 
dynamics and hence the transport properties of Eq. ([T]). 
The parameter xq is appropriately chosen such that the 
minima of V{xi^X2) are located at the integers; and S = 
1.614324, C = 0.0173, xo = 0.82, b = 0.1 and uj = 0.67 
are kept constant throughout. 

The 2D ratchet potential ([2]) is shown in Fig. [1] for four 
different values of the coupling strengths. The minima 
and maxima of the potential are marked with blue and 
red colors, respectively. It is to be noted that as k is 
increased, the heights of the potential V{xi^X2) move 
outward, opening up a valley along the diagonal, in which 
the two interacting ratchets may most likely share. 

Fig. [2] displays the behavior of a single trajectory anal- 
ysis. Here, on-off intermittency (001) [18] can clearly be 
observed in Fig. ^a). Fig. [2l[c) is the zoom of Fig. 
[2](b) that shows typical trajectories of the system. Dur- 
ing the earlier dynamics {t < 400), the particles advance 
in an asynchronous manner and then get captured in 
the ratchet potential well, where they achieve temporal 
synchrony - moving quasiperiodically as a single entity 
for a while, and suddenly escapes from the well and re- 
sume chaotic motion. The temporal synchrony state in 
which they move with the same displacement along the 
same path, is short-lived and then followed by unlocking. 
Thus, the repeated sequence of this motion gives rise to 
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FIG. 2: (Color online) Typical on-off intermittency (a) of 
the error state Ax{t) = xi{t) — X2(t) for k = 0.065, (b) of 
trajectories for k — 0.065 and k — 0.45 and (c) is the enlarged 
portion of (b) (300 < t < 800) for k = 0.065. 



001 synchronization. To quantify the 001 observed here 
we have computed the probability distribution of the lam- 
inar phases A(t) of Ax(t) and the average laminar lengths 
(/) of the trajectories xi,2(^) as function of e = a — ac at 
the critical bifurcation point ac = 0.0809474 as shown 
in Fig. [3l In Fig. [31(a), the collective dynamics for the 
coupled ratchets shows a —3/2 power law scaling, typical 
of on-off intermittency, while each subsystem dynamics 
exhibits type I intermittency similar to the single ratchet 
dynamics [19] with a —1/2 power law scaling as shown 
in Fig. [3^b). The above picture based on a single tra- 
jectory analysis is exact for single attr actor systems and 
may turn out to be misleading for irregular one, where 
periodic and chaotic attr actors could co-exist. 

Next, we explore the dynamics of the coupled ratch- 
ets ([1]) as is varied. Considering that single trajec- 
tory dynamics will not suffice for a highly chaotic sys- 
tem, all observables have to be averaged out over a large 
number of trajectories generated from the entire space 
[—1,1] X [—1,1] which is the unit cell of the resulting 
periodic structure. Here we make use of two important 
indicators, namely the error state 77 as a good measure of 
the synchronization and the current J as the transport 
quantifier. For a long time dynamics T, the error state 
for a given trajectory is given by: 



Vj 



1/2 



{X2^^^ - xi^^^)^ + (4^'^ - x[^^) dt, (3) 



with the full error r] = N ^"^j^i^Vj^ evaluated over the 
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FIG. 3: (Color online) Distribution of (a) laminar phases A(t) 
of Ax(t) satisfying a —3/2 power law scaling typical of on-off 
intermittency and (b) average laminar lengths with varying 
parameter a satisfying the scaling law (/) oc e~°'^ with e = 
ac — a and ac = 0.0809474 [xi{red) and X2{green)] showing 
type-I intermittency. Here, k = 0.065. 



total number N of trajectories. On the other hand, the 
current in a subsystem {i = 1,2) is defined as follows: 



Ji 



1 



M -UrN 



^ M N 

EE 



(z = l,2). (4) 



^=1 j=i 



where N is the total number of trajectories and ti of Xi{ti) 
a given observation time. This gives the average velocity, 
which is then further time-averaged over the number of 
observations M - Uc being an empirically obtained cut-off 
accounting for the transient effect, such that a converged 
current is obtained [Tl]. 

With the amplitude a = 0.0809472 of the driving, 
each independent system {k = 0) exhibits intermittent 
chaotic dynamics. Fig. [H^a) displays as a function of 
k in the range k G (0,1). Above the threshold value 
k > kth ^ 0.5, T] approaches zero, indicating a fully 
synchronized state. In Figs. [4l(b,c) we observe, in the 
same range of /c, the global dependence of the currents 
J_ = Ji — J2 and Ji, respectively, on the degree of syn- 
chronization. Fig. m^b) shows that prior to the fully 
synchronized state, J_ fluctuates around zero and when 
a full synchrony is achieved, J_ is identically zero. No- 
tice also that for 0.259 < k < 0.295, J_ is identically 
zero. The zero current here is not however due to syn- 
chrony, but rather to no transport, see Fig. IH^c). This 
figure also shows several spikes or jumps in Ji prior to the 
synchronized regime, in which (i) current-reversals show 
up for lower values of /c, (These spikes with or without 
current-reversals certainly have dramatic changes in the 
corresponding bifurcation diagrams [ll^), (h) a regime 
of zero-current occurs in the range (0.259 < k < 0.295) 
and (iii) full synchronization takes place and is accom- 
panied by current-reversal, during which the current is 
marginally zero and then negative for all k > 0.51. This 
remains to be further investigated, if this transition with 
synchrony is a characteristics feature of such systems. 
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FIG. 4: (Color online) Transition to full synchronization 
showing (a) the average error dynamics, 77; (b) the current, 
J- = Ji — J2 and (c) the current, Ji in the same range of k. 
The parameters are set as: a = 0.0809472, 6 = 0.1, xo = 0.82, 
and 00 = 0.67. 
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FIG. 5: (Color online) Current J vs the amplitude a for (a) 
strong couplings, 0.4 < /c < 1 and (b) weak coupling regime 
<k< 0.07 for fixed k. 
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Finally, special attention is markedly paid to control- 
ling current reversals as these may happen to be unde- 
sirable as far as transport is concerned. The case k = 
corresponds to no interactions for which current-reversals 
have been observed in a single ratchet model [1^, llH- 
The large coupling regime, where full synchronization 
is reached, corresponds identically to the current rever- 
sal observed in one single rachet model, (see Fig. [5](a)). 
However as k takes on smaller values, dramatic changes 
occur on current reversals leading to the rectification of 
the particles motions - for instance, current reversal ob- 
served in a large window of k (Fig. [5] (a)) is completely 
eliminated in Fig. E^b). Likewise current reversals found 
for smaller values of k are destroyed as k increases. Inter- 
estingly, there are few values of the coupling strength for 
which the system becomes totally reversals- free. This re- 
sult clearly demonstrates the importance of the coupling 
strength for the full control of transport over currents. 

To sum up, we have shown that the dynamics of a sin- 
gle particle ratchet can be significantly modified when 
coupled elastically to a second one. We have thus made 
use of the coupling strength to systematically rectify the 
particle motion. In particular, current reversals observed 
in a single ratchet can completely be annealed by appro- 



priately choosing the coupling strength. Interestingly a 
current-reversal-free regime has been also detected. The 
particle-particle interaction thus plays a significant role 
in the transport properties and acts as a control input 
that can be used to track current reversals. It thus 
implies, that the nature of particle-particle interactions 
could be used to control the transport properties of non- 
equilibrium dynamical systems. Moreover we found, that 
for larger values of the coupling strengths, the dynamics 
of the single particle ratchet predominates, leading to the 
fully synchronized state. 
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